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We study the 0{N) scalar field theory with quartic self-coupling in de Sitter space. When the 
field is light in units of the expansion rate, perturbative methods break down at very low momenta 
due to large infrared logarithmic terms. Using the nonperturbative large- A'^ limit, we compute the 
four-point vertex function in the deep infrared regime. The resummation of an infinite series of 
perturbative (bubble) diagrams leads to a modified power law which is analogous to the generation 
of an anomalous dimension in critical phenomena. We discuss in detail the role of high momentum 
(subhorizon) modes, including the issue of renormalization, and show that they influence the dynam- 
ics of infrared (superhorizon) modes only through a constant renormalization factor. This provides 
an explicit example of effective decoupling between high and low energy physics in an expanding 
space-time. 

PACS numbers: ll.lO.-z, 04.62. +v 

Keywords: Quantum field theory, de Sitter space, large-A teciiniques 



Ph. 



> 

m 

O 

cn 



X 



I. INTRODUCTION 

The study of quantum field dynamics on de Sitter 
space has received strong phenomenological motivations 
with the impressive observational success of the inflation 
paradigm and the need to compute quantum corrections 
to inflationary observables and with the observation of 
the recent acceleration of the Universe. Radiative cor- 
rections on de Sitter space have been addressed in a va- 
riety of field theories, mainly based on perturbative loop 
expansion [11-01 • the phenomenologically relevant case 
of a light field, with a mass small as compared to the 
Hubble parameter, loop diagrams typically exhibit large 
infrared (IR) logarithms, unless they are protected by 
some symmetries [8l-[l]|. In addition, in cases where the 
de Sitter symmetry is broken, e.g., by finite time ini- 
tial conditions, loop corrections give secular divergences 
which grow as powers of the number of e- folds [1, [ll| - [T5| . 
The understanding of such IR/secular divergences is of 
key importance for a fundamental description of quan- 
tum field dynamics in de Sitter space. For instance, it is 
unknown whether they signal an instability of de Sitter 
space against quantum fluctuations fl6l - l22j . In any case, 
such divergences pose a problem of principle and call for 
resummation. 

The situation is somewhat similar to the appearance 
of large logarithms in standard quantum field theory 
(QFT). These signals a breakdown of perturbation the- 
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ory and their resummation lead to nontrivial phenom- 
ena such as running couplings or anomalous scaling of 
field correlators [H, [^^l. IR and/or secular divergences 
are common features of perturbative approaches e.g. for 
bosonic theories at high ternperatures [25| or near a sec- 
ond order phase transition [2J] or in the case of nonequi- 
librium systems (26| . Various resummation techniques 
or nonperturbative approaches have been developed over 
the years to cope with such issues in fla t sp ace-time. 
These include, e.g., large- iV techniques [U [2^, the 
renormalization group [2J, [2^ , hard thermal loops [25| , 
or two-particle- irreducible (2PI) techniques [l^ls^l- Re- 
cently, some effort has been put in trying to extend these 
techniques to the case of IR/secular divergences in de 
Sitter space [ll,[3lH4l|. 

The large-iV limit in 0{N) scalar theories provides 
a simple nonperturbative approach which captures non- 
trivial IR physics in flat space-time [13, [H, iSjil ] . Re- 
cently this approach has been _applied in de Sitter space 
with interesting results [Slj, |33|. For instance, it de- 
scribes dynamical mass generation, first pointed out in 
[T^ . which screens perturbative IR divergences and ren- 
ders the theory well behaved. The generated mass ex- 
hibits a nonanalytic dependence in the coupling, which 
is typical of nonperturbative IR physics. Also it has been 
shown that strong IR fluctuations prohibit the possibility 
of a spontaneously broken phase [33l l45j , a phenomenon 
akin to what happens in two-dimensional flat space-time 
[3. 

The phenomena described here involve the resumma- 
tion of a particular class of IR divergences correspond- 
ing to purely local contributions to the self-energy. Such 
masslike contributions are specific to theories with quar- 
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tic self-interactions and are technically easy to deal with 
due to their local character. Less simple but more generic 
are the logarithmic IR divergencies arising from nonlocal 
perturbative diagrams. 

In the present paper we consider the four-point vertex 
function of an 0{N) scalar field theory. In the large- iV 
limit, the latter is given by an infinite series of nonlocal 
bubble loop diagrams, each of which bringing additional 
IR logarithms. Exploiting the physical momentum space 
representation (hereafter called p-representation) for de 
Sitter correlators [39, 113, 13, we show that this series 
can be exactly resummed. This induces a modification 
of the IR power law behavior in a similar way as the 
resummation of large logarithms in standard QFT, e.g., 
at a critical point, produces an anomalous scaling of field 
correlators. 

The nonlocal character of the physics described above 
leads to an nontrivial interplay between IR and ultra- 
violet (UV) momentum modes. In the main body of 
the paper, we simply disregard high-momentum modes. 
This allows us to exhibit the anomalous dimension phe- 
nomenon in arbitrary dimension. We show that the re- 
sults are essentially insensitive to the scale which sepa- 
rates IR and UV modes. A detailed analysis of the role 
of UV modes in four space-time dimension — including a 
discussion of renormalization — is presented in the Appen- 
dices together with some technical material needed for 
the main body of the paper. Remarkably, we find that 
UV modes do not alter the IR power law behaviors and 
merely give rise to a finite computable constant renor- 
malization factor. This demonstrates a form of effective 
decoupling between IR and UV physics in de Sitter. 

In Sec. In] we briefiy review the main aspects of the 
p-representation of de Sitter correlators. Then Sec. IIIII 
focuses on the 0{N) scalar theory in the large-TV limit 
for both the two- and the four-point vertex functions. We 
recall the solution of the corresponding equations for the 
two-point function, which describes self-consistent mass 
generation in Sec. IIVI The analysis of the four-point 
function and the resummation of nonlocal IR logarithms 
is performed in Sec. |Vl where we neglect UV modes. 
The important technical details needed for this analysis 
are presented in Appendices [AllCl Finally, Appendices 
iDllO present a detailed treatment of UV modes in four 
space-time dimensions, where the present theory is per- 
turbatively renormalizable. 



II. p-REPRESENTATION OF TWO-POINT 
CORRELATORS IN DE SITTER SPACE 

We consider a generic quantum scalar field ip{x) on the 
expanding Poincare patch of de Sitter space-time with 
D = d + 1 dimensions and expansion rate H = 1. The 
line element is given by 
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FIG. 1: The closed path C — C U in conformal time rj. 
The forward (upper) branch goes from — oo to 0~ and the 
backward (lower) branch C~ goes back from 0~ to — oo. 



in terms of the conformal time rj and comoving spatial 
coordinates X, with a(ry) = — l/rj = with t G R. 

Like general nonequilibrium quantum systems [26l . |49| - 
Is^ . quantum field theories in cosmological spaces can 
be conveniently formulated using a closed contour C in 
the time coordinate [H, [H, [s^] • The appropriate con- 
tour for conformal time is depicted in Fig. [T] The vari- 
ous components of 7i-point correlators are then described 
by means of time ordered products of field operators 
on the contour. For instance the two-point function 
G{x,x') = {Tc^{x)ip{x')) , where Tc denotes time or- 
dering on the contour C, encodes both the statistical 
and spectral correlators^ F{x,x') = ^{{ip{x),(p{x')}) and 
p{x,x') = i{[Lp{x),if{x')]): 

Gix.x') = F{x,x') ~ ^signc(a:'' - x'°)p(x,x'), (1) 

where the sign function is to be understood along the 
contour C. In the rest of this subsection, we consider the 
statistical function F . Everything holds equally for the 
spectral function p. 

De Sitter invariance ensures that the two-point func- 
tion F(x^x') only depends on the invariant distance 
z{x, x'). In the coordinate system ([1]), the latter reads 



z{x, x') 



+ T^'^ - (X - XQ 
27777' 



l\2 



(2) 



The dependence on comoving spatial variable is only 
through |X — X'l, which refiects the spatial homogeneity 
and isotropy of de Sitter space-time in this coordinate 
system. 

It proves convenient to introduce conformally rescaled 
quantities, such as the rescaled field 4>{x) — {r/)ip{x) 
and its correlators. One has, for instance, 

F{x, x') - [a(7;)a(7;')]"'^ ^,(77, 77', |X - X'|), (3) 

where i^c(f7,?7MX- X'l) = i({0(a;), ^(a;')}). Introducing 
spatial comoving momentum variables, one defines 



FM.K)^ d''Se'^'^-^F,{r^,r,',\S\ 



(4) 



Spatial homogeneity implies the conservation of the co- 
moving momentum. 



a^iv) {-drf+dX.^) 



1 We use {A, B} = AB + BA and [A, B] = AB - BA. 
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FIG. 2: The closed path C = C+ U C~ in the momentum 
variable p — —Ktj. The upper branch C'^ goes from +00 to 
0^ and the lower branch C~ goes back from 0"'' to +cx3. 



Specializing to de Sitter, an extra symmetry imp lies 
that two-point correlators follow the scaling law |39l, llsj 



(5) 



where p — —Krj and p' — Kr]' are the physical momenta 
at times 77 and 77' respectively. Eq. ([S]) provides the p- 
representation of the two-point correlator. 

As mentioned previously Eq. ([5]) applies to the spec- 
tral function p as well. In fact, as in ([Ij, one can combine 
the statistical and spectral two-point functions in a sin- 
gle propagator defined on a closed contour C in physical 
momentum 1391: 



G{p,p') = F{p,p') - -aign^ip- p')p{p,p'). (6) 

The sign function is to be understood on the contour 
depicted in Fig.O We also note the symmetry properties 
F{p,p') = F{p',p) and p{p,p') = -p{p',p). 



III. 0{N) SCALAR FIELD THEORY IN THE 
LARGE- iV LIMIT 



We consider an 0{N) scalar theory, defined by the 
classical action 



fa- ^ifa^a)'^ } , (7) 



where a = 1, . . . , N and a summation over repeated in- 
dices is understood. Here, 



□ ; 



1 



is the covariant Laplace operator and 



MS 



(8) 



(9) 



is the effective square mass with m the tree-level mass 
and ^ the coupling to the Ricci scalar R = d{d +1). 
Finally, 



= / d^x^/~g{x) = J dx" I d'^xy/-g{x), (10) 



where the a;*^ integral runs along the contour C. 




FIG. 3: The self-energy T,{x,x') in the limit N ^ oo. The 
internal line in the diagram is given by the full propaga- 
tor G{x, x), hence the nonperturbative nature of the large- 
limit. This describes self-consistent mass generation. 



A. The large-iV limit 

The 1/N expansion provides a controlled expansion 
scheme which gives access to nonperturbative aspects of 
the theory. It can be formulated in cosmological spaces 
[39! . [ssj using standard functional techniques, provided 
one employs appropriate covariant quantities such as the 
functional derivative. 
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Sip{x) ^/~g{x) 5(p{xy 



(11) 



the integration measure (jlOl) and the corresponding co- 
variant Dirac distribution on the contour. 



(12) 



defined such that J S'-^^^x, z)f{z) = f{x) for any func- 
tion / on the contour. Let us recall the basic equations 
describing the limit N ^ 00, see [H, HJ. It has been 
shown in [33[ that the theory does not admit a de Sit- 
ter invariant state with spontaneous symmetry breaking. 
We thus consider the symmetric phase, where {(fa) — 0, 
Gab = SabG and Sab — (5afcS. The covariant inverse prop- 
agator is defined as 

G-\x,z)G{z,x') ^ 6'-^^x,x'). (13) 

It can be written in terms of the covariant self-energy as 
G-^{x,x') ^Gq\x,x')-Y.{x,x'), (14) 

where 

tG^\x,x') = (□,-m2s)<S(^)(x,x'). (15) 

In the large- iV limit, the self-energy is given by the tad- 
pole diagram of Fig. [31 where the line represents the 
propagator G itself. It thus corresponds to a local, 
masslike contribution; 



where 



i:ix,x') = -iaS'-^\x,x'), 



o- = -G{x,x). 



(16) 



(17) 
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FIG. 4: The infinite series of bubble diagrams contributing 
to tlie function I{x,x'). The black dots denote interaction 
vertices whereas the crosses are the end points of the func- 
tion. The elementary one-loop bubble is given by the function 
Tl{x,x') (X G^{x,x'). Each addition of a new bubble leads to 
large IR logarithms in the momentum representation, the in- 
finite series of which resums to a modified power law. 



De Sitter symmetry guarantees that the latter is a con- 
stant. The propagator G thus satisfies the equation 

(D,- M^)Gix,x') ^iS^^\x,x'), (18) 

with a self-consistent square mass 

= mis + ^- (19) 

The four-point vertex function reads, writing Xi = 
xi, . . . ,xa for brief, 

ritL(^0 = 5ab5cd5''°\xi,X2)5''°\x3„X4)iD{xi,X3,) 

+ perm., (20) 

where "perm." stands for the two cyclic permutations 
needed to make r|,tcd(^») symmetric. Here, the function 

iD{x,x) ^ 5^^\x,x')+il{x,x) (21) 

is the two-point correlator of the composite field x — 
\ipaipa/6N: D{x,x') = {Tcx{x)x{x')). The local contri- 
bution gives rise to the classical vertex in ([20)1 and the 
nonlocal one actually corresponds to the infinite sum of 
bubble diagrams represented in Fig. 2] This infinite sum 
can be cast into the following integral equation: 

/(x, x') = n(a;, x')+i / n(a;, z)/(z, a;'), (22) 



with 



Ji{x,x') = ~-G'^{x,x') 



(23) 



the elementary one-loop bubble diagram. The four-point 
vertex is thus completely expressed in terms of a nontriv- 
ial two-point function. 



B. p-representation 

The p-representation plays an important role in the 
analytical solution presented below. We briefly review 
the relevant material for our present purposes and refer 
the reader to Ref. ISSl for more details. 



Eq. (|18l) is equivalent to the free-field-like equation for 
the statistical function F defined in Eqs. (IS])-® 
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where 



M2. 



(24) 



(25) 



The spectral function p{p,p') satisfies the same equa- 
tion. The "initial" data is to be specified a.t p = 
p' — >■ oo. Commutation relations imply p{p,p')\p=p ~ 
dpdp'p{P^P')\p=P' = and dpp{p,p')\p^p> = -1 for the 
spectral function. The statistical function contains the 
information about the (quantum) state of the system. 
The only choice compatible with de Sitter symmetry and 
with the requirement of renormalizability is the Bunch- 
Davies vacuum This corresponds to the choice 



F{p,p') 
dpF{p,p') 
dpdp,F{p,p') 



1 

2' 
0, 
1 

2' 



(26) 



With this choice, the solutions of Eq. ([M)) for F and p 
read 



F{p.p') = -^Kg{H,{p)H:{p')} 



(27) 



/^lm{H,{p)H;{p')}, (28) 



where Hi, (x) is the Hankel function of the first kind |5(: 
Finally, the self-energy pT|) reads 



A 
6 

d''q 



F{q,q) 



Att 
24 



\HA<1)\' 



(29) 



where /q = / 
solution of the gap equation ([T9l) with (|29|) and (|25|) . For 
instance, in the case of vanishing or negative tree-level 
(classical) mass m^g < 0, where perturbation theory is 



The self-consistent mass M is the 



ill-defined, Eq. (|T9l) describes the dynamical generation 
of a strictly positive square mass which cures the IR di- 
vergences of perturbation theory. 

Let us now consider the four-point vertex function (|20p . 
Introducing the appropriate confornial factors 



r^^' (x- 



[a(?7i)-- •0(774)] 



T.(4) 



cabcd 



ix^) (30) 



and going to spatial comoving momentum space, we have, 
extracting a factor {2tt)'^ S^'''> (j2t=i K 



;rf(7?i>Ki) = [0(771) • • •0(774)]" 



: \^SabScd Sc ivi - m)Sc im - rii)^Dc{r]i , 773, iv:i2) -l-perm.| , 

(31) 
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where i^ij = |Ki+Kj|. Here, Dc is the comoving momen- 
tum representation of the function D, Eq. (PT|) . defined 
as 



Dix,x')^[a{7j)a{7j')] - r/MX - X'|) (32) 



and 



D,ir^,rj',K)^ / S e-^"^-"- D4rj,Tj' ,\S\) 



(33) 



We define the comoving momentum representations Ic 
and He of the functions / and 11, Eqs. (f22|) and (|23|) . in 
the same way as for Dc- Equation (|2ip then reads 

z^e(77,r?',/0 = - ^ [6c{v-v')+^Uv,v',K)] . (34) 

As aheady emphasized, the function Ic encodes the non- 
trivial (loop) contribution to the four-point vertex in the 
large- limit. It corresponds to the infinite sum of bub- 
ble diagrams depicted in Fig. |4] with the elementary one- 
loop bubble He- The functions Dc and He and Ic admit 
the following p-representation: 



DciT],v',K)^KD{p,p') 



(35) 



and similarly for He and /c, where p = —Kr] and 
p' — —Kr]'. Introducing the Dirac distribution on the 
momentum contour^ C 



one has 



6c{v - v') = -K6Ap - p'), 



A 



iD{p,p') = -— ^-S^(p-p') + il{p,p') 



The one-loop bubble 11 is given by 

n(p,p') = 



A^^^,^£^ f G{qp,qp') G{rp,rp') 
6 Jn q r ' 



(36) 
(37) 

(38) 



where r = |e + q|, with e an arbitrary unit vector, and 
the function / solves the following integral equation on 
the contour C: 

iip,p')=Il{p,p')-i f dsliip,s)iis,p'). (39) 
Jc 

Finally, we write the above equations explicitly in 
terms of their statistical and spectral components on the 
contour, defined as'^ 



n(p,p') = flF{p,p') - ^sigiicip - p')Up{p,p'), (40) 



with Uf{P:P') = Uf{p',p) and np{p,p') = -Up{p',p) 
and similarly for /. One has 



and 



'^iPP') 



— \F{qP,qp')F{rp,rp') 
qr L 

^/5(gp, qp')p{rp, rp') | (41) 



\{PP') 



F{rp, rp') p{rp, rp') 



qr 



(42) 



Writing explicitly the integrals on the contour C in p9p , 
one obtains the following equations 



ip(p,p')^fip(p,p')+ dsnp{p,s)ip{s,p') (43) 



and 



iF{p,p') - flH{p,p') + / dsUp{p,s)iFis,p'), (44) 
Jp 

where we introduced the auxiliary function 

/•OO 

nH{p,p')^nF{p,p')- dsiiF{p,s)ip{s,p'). (45) 



For later use, we note that, using (|43l) . Eq. (|44| can be 
solved in terms of the function Ip as, see Appendix [XI 

iF{p,p') = Uh{p,p') + / dsip{p,s)flH{s,p'). (46) 
Jp 

In the next two sections, we solve the above equations 
in the deep IR limit, that is for p,p' <C I in the case of 
light fields, <^ 1. We adopt the following general 
strategy: we first assume that the IR physics is domi- 
nated by IR modes and we cut all momentum integrals 
(either in momentum loops or in "time" integrals) at a 
momentum scale /i < 1. We shall see that we can find a 
closed solution of the equations in that case and we ver- 
ify a posteriori that the solution weakly depends on the 
scale /i. Of course this is not a proof of our assumption, 
merely a consistency check. We defer to the Appendices 
a more complete analysis taking into account the role of 
high momentum modes. There we show that our results 
only get slightly modified by a constant renormalization 
factor. Similarly, we neglect renormalization aspects in 
the main course of the text and defer their discussion to 
Appendix [G] 



^ The latter is defined sucii tliat dp' S^(p — p')f{p') = f{p) for 
any function / on the contour C. The minus sign in I I36I I copes 
for the "wrong" orientation of the momentum contour. Note, in 
particular, that dp = —Kj^ drj [39ll . 

3 This decomposition is modified by renormalization. However, the 
final equations to be discussed here are unchanged, see Appendix 

m 



IV. TWO-POINT FUNCTION: MASS 
GENERATION 

The solution of the gap equation (|T9)) has been dis- 
cussed in [31I - I33I . Here, we briefly review the main re- 
sults and introduce some necessary material for the next 
section, where we compute the four-point function. 
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In the case of light fields, with M < d/2, the index v 
is real and < < d/2. It is customary to introduce 
the parameter e — d/2 ~ v. For small mass, one has 
e « ]\/P /d. Following the strategy described above, we 
only need the statistical and spectral correlators F and p 
in the IR regime, i.e., for p,p' ^ 1. The latter are given 
by, see Appendix iBl 



{pp'Y' 



where = [2T(i^)]V47r and 



2v 



It proves useful to rewrite the latter expression as 



In^ 
p' 



with 



(47) 



(48) 



(49) 



(50) 



Eqs. dHj) and g!]) by their IR expressions (|471)-(150l). The 
loop integrals are performed in Appendix [UJ We obtain 



with K = I — 2e (notice that n ^ v — e < v) and 



tvfip^p') = ^Vti\n^ 
vPP \ P 



where we defined the odd function 
with Vu{x) given in (j50p . and where 



6e {2nY 



2a and ttf 



-T^pFy. 



(54) 



(55) 



(56) 



(57) 



Let us first consider Eq. for Ip . We seek a solution 
of the form 



Ifip.p') 



/pi/ 



(58) 



Restricting the momentum integration in ([29)) to the 
IR part |q| < /i, one gets 



where 



A 

"^6 



FiR((7,g) XF, 



12e (27r) 



(51) 



up to corrections of relative order eln/i ^ 1. Here fl^ = 
27r^/^/r(d/2). Using e w M^/d in the small mass case, 
it is easy to solve the gap equation (fTO)) . For instance, in 
the case of vanishing tree-level mass, m^g — 0, one gets 



IXF, dQd 
12 {2nY 



so that 



(T 

d 



a{d/2) 
127r 



(52) 



(53) 



where a{x) = ^J'iT{x) / xtt''^^ is chosen here so that 
a{i/2) = 1 for c? = 3. The self-consistent square mass 
(|52p is nonanalytic in the coupling, reflecting its nonper- 
turbative IR character. 



V. FOUR-POINT FUNCTION: RESUMMING 
INFRARED LOGARITHMS 

The first step in the calculation of the four-point vertex 
is the evaluation of the bubble n, Eqs. (gT]) and (|42l) . 
Restricting the loop integral to the contribution from IR 
modes, i.e. such that qp, qp' < /i, which also implies 
rp, rp' < /z, one can replace the propagators F and p in 



lix). 



(59) 



with T{x) an unknown odd function, as implied by the 
antisymmetry of Ip{p,p') in the exchange p ^ p' . The 
exp(— e|a;|) terms cancel out in Eq. (j43|) and we are left 
with the following integral equation: 

I{x) = r,{x) -TTp f dyV,{x - y)I{y). (60) 
Jo 

As a first attempt to solve this equation, one may try 
an expansion in the parameter TTp which counts the num- 
ber of bubbles. This is akin to perturbation theory. At 
first order, one has 

lix) = r,{x)-Trp f dyVu{x - y)V.{v) + O {^D 



tanh(i/a;) 



(61) 



For \x\ ^ 1, that is for large logarithms or, equiva- 
lently, large cosmological time separations, |ln(p/p')| = 
|t - t'l > 1, we get 



I{x)^rAx) {l~^\x\+0{7r^ 



(62) 



Thus, we see that the true expansion parameter is TTp\x\ — 
TTpl ln(p/p')| and that the perturbative solution breaks 
down for | ln(p/p')| = |t — t'| > I/tt^. The argument ex- 
tends to higher orders and one can check that, in that 
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case, one has to resum the infinite series of bubble dia- 
grams. 

Fortunately, the exact solution of Eq. (j60| is easily 
obtained. It reads 



I{x) 



sinh fx 



where (note that v < v) 



We thus have 



/pp' 



P 



(63) 



(64) 



(65) 



We see that resumming the infinite series of bubble 
diagrams — each producing higher powers of ■Kp \ lnp/p'| — 
through the integral equation (l43l) changes the exponent 
V characterizing the IR power law behavior of the single 
bubble function Hp, see Eq. (j55p . into the new exponent 

(IMI). 

Let us now consider the case of Ip, Eq. ([33]), where, 
following the general strategy described previously, we 
replace the upper bound of momentum integrals in (j44l) - 
(|46p by /i < 1. The contribution of higher momenta is 
analyzed in Appendix |F1 We first note that the F compo- 
nent of the symmetric bubble function 11 can be written 
as the following factorized form 



= In^U In^ 



(66) 



where we artificially introduced a dependence on /i for 
later convenience. Here, 



A{x) 



(67) 



The factorization property (|66l) together with the fact 
that Eq. (HSI) for n^(p,p') only involves the function 
n^ (p, s) with fixed first argument imply that the auxil- 
iary function takes the asymmetric factorized form. 



"^^fln^ 



A In 



/i 



(68) 



with the function A given in (1571) above and the function 
A is to be determined. Using (j45p and ((58)) . we find (here 
x = lnp//x < 0) 

A{x) = A{x) +Trp [ dyl,{x - y) A{v), (69) 



Finally, using the solution (|46|) of Eq. (|44)) , where the 

function n^(s,p') always appeasr with the same second 

argument as /J?(p,p')) "we conclude that the latter also 
factorizes: 



A\\n^ 



A[\n^ 



(71) 



with the function A given in Eq. (|70p above. The fac- 
torization property ((7T|) of /|P follows directly from the 
factorization of H^, Eq. ((66|) . As described in Appendix 
[Cl the latter follows from the corresponding factorization 
property (|T7)) of the statistical correlator /jr and is thus 
reminiscent of the fact that IR de Sitter fluctuations can 
be described as a classical stochastic ensemble of growing 
modes^ [57h59i] . 

Let us discuss in more detail the above solution, 
Eq. ((7T|) . Recalling that v < one finds, for —x ^ 1, 



Aix) 



2v{v — v) fi'^ 



(72) 



where we neglected terms of relative order 0{TTp) in the 
second equality. Here, we introduced the new exponent 

R = V — e. (73) 

Again, we find that the infinite summation of bubble di- 
agrams changes the IR power law of the single bubble 



A In 



1 



A In 



1 



(74) 



and K — K can be regarded as an anomalous dimension. 
We thus have 



{pp') 



(75) 



We emphasize that the dependence of the above results 
on the scale /i < 1 is quite mild oc /i"""^ « 1 — 2^ In/i, in 
the case iTp ^ 1. This suggests that the neglect of high 
momentum modes in the present analysis is a consistent 
approximation. 

As in the case of Ip, it is instructive to compare the 
exact solution (|72p to the perturbative expansion in pow- 
ers of TTp. Using ([SO]) . Eq. ([S^ reads, at lowest nontrivial 
order, 

A{x) = A{x)+7:p / dyreix-y)A{y) + 0{Trl) 



^(.)(1 + ^ 



^2vx 1 



2v 2v 



(76) 



where we used Ie{x) 



A{x) = ^A{x) ' 
2v 



4(-x). We get 



(70) 



* The fact that the connected two-point function factorizes is also 
to be put in parallel with the phenomenon of undulation in (ana- 
log) black hole physics [6Cj . In the present case, the undulation 
occurs in physical momentum space. 
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For IR momenta, —x = \n(fi/p) ^ 1, we get 



A{x) « A{x) [l + ^x + O [ttD 



(77) 



We see that the true expansion parameter is, in fact, 
TTphifx/p and that, for momenta ln(/i/p) > l/'^'p, one 
needs to resum the IR logarithms generated by the series 
of bubble diagrams. This results in the modification of 
the IR exponent k k. 



VI. DISCUSSION 

In the previous sections, we have mainly considered 
the case of a small, or vanishing, tree-level square mass 



MS 



strong IR divergences. As previously emphasized, one 
has, in that case, e oc a/A and thus TTp cx VX. The actual 
expansion parameter is thus \/A. The mass generation 
phenomenon discussed in Sec. IIVI can be viewed as a 
sort of anomalous dimension in the IR [36| . Indeed the IR 
power law behavior of the free propagator of a massless 
field, characterized by the exponent d/2, is changed to an 
exponent v d/2 — e which regulates the IR divergences 
of the perturbative series. 

Similarly, the resummation of nonlocal IR logarithms 
in the two-point function of the composite operator ip"^ 
or, equivalently, in the four-point vertex function of the 
field (f leads to a modified IR behavior: the power law 
behaviors of the one-loop bubble contributions Hp and 
n^, characterized by the exponents v ~ d/2 — e and 
K « d/2—2e, are changed to « d/2~2>e and R « d/2~Ae 
respectively. 

We point out that, inversely, the dynamically gener- 
ated IR power law behavior of the function Z), for large 
(super-Hubble) momentum/time separations |lnp/p'| = 
\t — t'\ > 1, can be viewed as the generation of an effective 
mass for the composite field ip'^. Using, see (l73l) . 



^I^^l sinh; 



sinh 



tiX ^ 



for 



>1, 



(78) 



as well as the relation (j57p between -Kp and TTp, one has, 
for I \np/p'\ > 1, 



(79) 



pp'Df{p,p') « ^^^V-.(\n^)j . (80) 

We have assumed e <C 1 and TTp <C 1. In this case, one 
can identify R k D k v in the constant factors on the 
right-hand side. Using Eq. , we see that the function 
pp' D^^{p,p') behaves as the two-point function of a free 
massive field with square mass = — k^. For 
vanishing tree-level mass, we get M ^ 2M. 

There is another interesting situation to be mentioned, 
namely the case of a negative tree-level square mass 
< 0. This is the case where the classical potential 



shows spontaneous symmetry breaking. In that case, per- 
turbation theory is again ill-defined because of the pres- 
ence of massless Goldstone modes. It has been shown 
in '33'] that, as in the previous case, nonperturbative ef- 
fects of self-interactions cure the theory by generating a 
positive mass square restoring the symmetry. The self- 
consistent mass is approximately given by 



MS 



from which it follows that 



4 12 {2'kY 12|m2g| (27r)'i' 

(81) 

(82) 



— )■ 0. In this case, perturbation theory suffers from and thus 



TTp = 2\m^^\ 



(83) 



m 



dS 



This is equal to the radial curvature of the tree-level po- 
tential at its minimum. We see that, in that case, e cx A is 
parametrically smaller than in the case of vanishing tree- 
level mass, leading to stronger IR enhancement and to a 
truly nonperturbative value of the parameter TTp govern- 
ing the bubble resummation. At weak coupling, we get 
M K, ^Jtt p ^ M . The composite field ip^ gets a nonper- 
turbatively large mass. 



VII. CONCLUSIONS 

Resumming perturbative IR divergences of QFT on de 
Sitter space is a timely issue of modern cosmology. The 
present study provides an exact solution of a nonpertur- 
bative limit where large IR logarithms can be explicitly 
resummed. This gives rise to modified power laws for 
the correlators in the IR, akin to anomalous dimensions 
in the context of critical phenomena. 

A detailed analysis of the role of UV modes reveals 
that they only affect the IR results by a multiplicative 
constant. This is a form of effective decoupling between 
UV and IR physics in de Sitter space. That such a de- 
coupling occurs in an expanding background is a highly 
nontrivial result. 

The results obtained here in the large- TV limit are new 
exact results in the context of QFT on de Sitter. Since 
the early work of Starobinsky and Yokoyama [l2|, non- 
perturbative calculations have been rather scarce, owing 
to the difficulty of extending to de Sitter space standard 
flat space-time resummation or nonperturbative tech- 
niques. The situation has improved in recent years with, 
in particular, the use of the dynamical renormalization 
group approach [15i], explicit QFT calculations in vari- 
ous nonperturbative approximation schemes [31. 33-36], 
all orders perturbative results [2l|, l2p. [61I , a nd resumma- 
tion of Schwinger- Dyson equations j32l.l40l liTI , l62j . Also 
worth mentioning is the interesting proposal for a reorga- 
nized perturbative expansion in the context of Euclidean 
de Sitter space proposed in [STllssj. 
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We believe the methods proposed here and in [39| pro- 
vide a powerful tool to pursue these studies. In particu- 
lar, they are well suited to the study of the IR properties 
of de Sitter correlators and are in principle applicable to 
other approximation schemes. Examples presently under 
study include corrections to the two-point function 
in the 0{N) model and approximations based on two- 
particle-irreducible techniques [63]. 
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Appendix A: Proof of Eq. (|46|) 

Consider the integral on the right-hand side of Eq. p6 
and use Eq. (HH) to eliminate YIh- 



dsIp{p,s)UH{s,p') ^ / dsIp{p,s)lF{s,p') 
p Jp 

OO /"CO 

ds / ds' ip{p,s)i].p{s,s')iF{s',p'). (Al) 

p J s 



Writing the double integral on the right-hand side as 

(A2) 

and using Eq. ([43)) . the second line of Eq. (|A1|) rewrites 
as 



OO nOO 



OO t^S 



ds I ds — I ds I ds 

Jp J p 



OO 



ds' / dsIp{p,s)np{s,s')lF{s\p') 



p •'p 

ds' |np(p, s') - Ip{p, ,s')} lF{s'y). (A3) 



We thus get 



dsip{p,s)TlH{s,p')= / dsilp{p,s)iF{s,p'), (A4) 
p Jp 



hence Eq. (|46 



Appendix B: Two-point correlators 

In this series of appendices, we discuss various tech- 
nical aspects of the present work. Appendices [B] and [C] 
present material directly used in the main body of the pa- 
per, in particular the calculation of the one-loop bubble 
functions and li^^, Eqs. (HDl-dnni). The remaining 
Appendices [DHG] provide a detailed discussion of the role 
of high momenta (Appendix IF|) . including the question of 
renormalization (Appendix[G|, in the case D = d+1 = 4, 



where the theory (O is perturbatively renormalizable. 
This requires the evaluation of the functions II™'p and 

n^^, defined below and discussed in Appendices ID] and 
IeI respectively. 

In the present section, we recall the form of the statis- 
tical and spectral two-point functions and their approx- 
imate forms in the IR and UV limits. As in the main 
text, we separate the two regimes with an arbitrary scale 
/i ~ 1. For a generic two-point function we in- 

troduce the notations 

V'^'ipJ) = Vip<f,,p'<^,), (Bl) 
I?u^(p,p') = P(p>/i,p'>M), (B2) 

for the cases where both momenta are IR or UV and 



2^mix 



{p,p') - Vip<fi,p' >fi), (B3) 

{p,p') = v{p>^Ji,p'<^i), (B4) 



for the cases where one momentum in IR and the other is 
UV. Notice that for (anti) symmetric functions 'D{p,p') — 
±V{p',p), one has = ±X'^™(p',p). 

The statistical and spectral two-point functions of a 
free scalar field of mass M in the Bunch-Davies vacuum 
are given by 

F{p,p') = ^VpP^e{H4p)H:{p')}, 



/^lm{H,{p)H:{p')}, (B5) 



with ly = d? 1^ — Iv'P. In the next sections we shall 
use the propagator ((6]) with both p and p' on the up- 
per branch of the contour C, where s\g\ig{p — p') = 
— sign(p — p') and thus 

G{p,p') = F{p,p') + -sign{p - p')p{p,p') for p,p' e C+. 

We recall the symmetry relations G{p,p') = G(p',p), 
F{p,p') = F{p',p) and p{p,p') = -p{p',p). 

We are concerned here with the small mass case M < 1 
for which the index is real. The leading low and high 
momentum behaviors of the Hankel function are given 

by m 

H,{p)^^(-) = -z,/^/„(p), for 
tiT \p J y np 

(B7) 

where the second relation defines the function and 
T 



hAp) 



, for p:^ 1, 



(B8) 



where ipu — ^{i^ + ^). 

The IR behavior of the statistical function is readily 
obtained as 



{pp'r^ 



(B9) 
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where F^, = (v)]'^ / An . The factorization property 
(|B9I) expresses the (approximate classical nature 

and the phase coherence of IR fluctuations in de Sitter 
(57I [59I . The leading IR behavior of the spectral function 
necessitates a more precise calculation. It can be written 
as 



Pir{p,p') = 



(BIO) 



where we introduced the function Vi, of (|50p . The mixed 
components are obtained from (jB7p . 



Fmix{p,p') = -fv{p) 

/5mix(p,p') = 2^(p) 

or, equivalently, 



sin(y - ip^) 

COs{p' - (fy) 

71 ' 



(Bll) 
(B12) 



Gmi^{p,p') = fu{p)- 



for p,p' e C+. (B13) 



Finally, the UV components have the following 
Minkowski-like expressions (with physical momentum in 
place of time) 



F\Jy{PiP') = ^cos(p-p') 
Pvv{p,p') = -sin(p-p') 
or, equivalently. 



(B14) 
(B15) 



Guv(p,p') = 



for p,p'eC+. (B16) 



In the deep UV regime p,p' ^ 1 and for subhorizon time 
separation |lnp/p'| — \t ~ t'\ ^ 1, one has p — p' k, 
—p{t — t') (recall that p'e* — pe*") and one recovers the 
(massless) Minkowski vacuum correlator 



Guv(p,p') e-*Pl*-*'l 



2p 



(B17) 



Appendix C: Calculation of ff^{p,p') 

Here, we compute the leading IR behavior of the one- 
loop bubble function Il{p,p'), defined in ((38|) . It is suf- 
ficient to consider the case where both p and p' lie on 
the upper branch of the contour and p < p' , such 
that sign^(p — p') — 1. We separate the momentum 
integral in three regions: low momentum modes, such 
that qp < qp' < /i; intermediate momentum modes 
with qp p, qp'', and high momentum modes with 

p^qp < qp'- 



1. IR contribution 

The contribution from low momentum, IR modes gives 
the dominant contribution. One needs to consider the F 
and p components, Eqs. (j4T|) and (|42|) separately. Here, 
the condition qp < qp' < ^ implies that rp < rp' < 
fi and both propagators under the loop integral can be 
approximated by their IR behavior, Eqs. (IB9p and (|B10p . 
In that case, due to the strong IR enhancement of the 
statistical propagator as compared to the spectral one, 
one can safely neglect the (quantum) contribution pp as 
compared to the (classical) one^ FF in lip. We get 



vviR/ ^/ f Pm{qP,qp)Fm{rp,rp') 
^fVP^p) ~ -Typp ) ' / 

1 



(CI) 



for the F component and 



P(p,p') « --{pp')^ ^ Fin{qp,qp')piR{rp,rp') 
^ 3 



qr 



1 



4(w')^^iR(FP')m(Fp') / (C2) 



for the p component, where we recall that r = |e -I- q| 
with e an arbitrary unit vector. Since the momentum 
integral in the second line of Eq. (jCip converges rapidly 
at large momentum {2v w d) we can safely ignore the 
upper bound. One cannot do so in the momentum in- 
tegral in the second line of Eq. (jC2p since the latter is 
sensitive to the upper bound. This results in a nontrivial 
p' /p dependence, as discussed below. 

The integral in Eq. (jCip is easily evaluated by using 
the technique of Feynman parameters. We write 



1 



T{2i 



(qr) 



2v 



dx 



[x{l-x)Y 



r2(i^)7o [xq^ + {1 - xy^] 



2v 



(C3) 



&wAxq'^ + {l — x)r'^ = [q+{l — x)e]'^+x{l~x). Performing 
the shift q + (1 — x)e — )• q under the q integration, one 
can perform all the angular integrations in d-dimensional 
spherical coordinates, to get 



i_ _ jT^iw f' r, q^-'jxji ~ x)r-' 
(2^)'^r2(^)y/^y„ [,2+^(i_:,)]2''^ 



{qr) 



2u 



(C4) 

where fid — 2i:'^/'^ /T{d/2). Finally, performing the 
change of variable u — q^ /x{l — x), the u and x inte- 
gration can be trivially performed. We obtain, recalling 



This is a typical classical (statistical) field approximation 0, IgsI . 
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v = d/2 



1 



_ r(f)r(f-2e) T^s) 

1 



(27r)'' £ 

The leading IR behavior is then 



l + 0{e) 



(C5) 



(C6) 



with n— ^ — 2e — V — e<v and 



TTp = 



6e {2TTy 



(C7) 



As aheady emphasized below Eq. (1551) . the function 
assumes a factorized form which is reminiscent of the 
classical coherent nature of IR de Sitter fluctuations. For 
later use we rewrite fl^ as 



iif{py)^A{p)A{p'), 



(C8) 



where the function A is related to the one introduced in 
([57| as 



A{p)^'-^=^ 



""^.^fln^ 



VP V 
The integral in Eq. (IC2[) is easily performed: 



(C9) 



ql<4 9 



which leads to 



»d 1 /m\ 
(27r)'' 2£ 



2e 



{p 



{2t:Y 2e 



1 + Oie) , (CIO) 



2. Mixed and UV contributions 

We now show that the intermediate and high momen- 
tum contributions to 11^^ are subleading in the limit 
<C 1. The contribution from intermediate momenta 
9P ^ M ^ QP' is such that q > jj,/p' ^ 1, which implies 
that r ~ g ^ 1 and thus rp < 1 < rp' . It follows that 
both propagators in Eq. ([55]) can be replaced by their 
"mixed" expression ()B13j) . We thus have to evaluate the 
integral 



<|q|<f 
Fir{p,p) 



^mix 

{qp,qp')Gn,i^{rp,rp') 
qr 



Qd Finip,p) ^g-2.(,p'-^,/,)^ (C14) 



(27r 



where we can safely replace v — > d/2 in the last line. The 
remaining integral is bounded: 



^p-2i(9p'-V.) 



< 



p dq 



P -P 



(C15) 



and one easily checks that the corresponding contribu- 
tions to either = Refl™ or nj," = -2ilmni^ are 
negligible as compared to those obtained in the previous 
subsection (notice also that there is no IR enhancement 
factor 1/e). 

Finally, the high momentum contribution can be esti- 
mated along similar lines. In Z) > 4 the integral diverges 
for p' ^ p and needs to be regulated. Here, we regu- 
late the theory by imposing a sharp cutoff A on physical 
momenta, i.e. by effectively replacing the propagators 
G(p,p') ^ G{p,p')6{K - p)e{A - p'), see Appendix [Gl 
We thus have fi < qp < qp' < A and, therefore, g 3> 1, so 
that r « g ^ 1 and rp' > rp > ^. Both propagators in 
([55)) can thus be replaced by their UV expression (|B16p . 
The corresponding contribution reads 



tifiP,p') 



where 



/pp' \p' 



(P,\V[\n^] for p<p', 



6e (27r)^' 



(Cll) 



(C12) 



Using the antisymmetry of the function Hp (p, p') , we have 
then 



ti^^iP,p')-^('-)vLlL 
Vpp \pJ \ p 



for p > p' . 

(C13) 



From this we get Eq. (|5 



Gvviqp, qp')GuY{rp, rp') 



<|q|< 



1 

4 Ja 

P 



qr 

J{q+r)[p-p'] 



<|q|<i 



qr 



dq 



4{2nY Jt. q^-'i 
Again an upper bound is given by 



o2ig(p-p') 



(C16) 



£ q' 



dq 



< 



li q' 



dq 

3-d 



(^) -(f) 



d-2 



(C17) 
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which is easily shown to give neghgible contributions to 
both and 11^^ in the low momentum limit p,p' ^ 1, 
at fixed cutoff A. 

We end this subsection with a more precise evaluation 
of the UV contribution (jC16p . useful for our discussion 
of renormalization in Appendix [Gl One has, for p ^ 1 
and X ^ p — p' finite, 



p dq 



hl q 



3-d 



3-d 



2ix 



3-d 2iiix 

p \ e p 



2ix 
(C18) 

We see that, as expected, that the one-loop bubble — 
which we recall is a one-loop contribution to the four- 
point vertex function — is finite in the limit A — >■ cxd for 
D = d + 1 < 3, where the present theory is super- 
renormalizable, and has no coupling divergence. The case 
D = 4, where the theory is (perturbatively) renormaliz- 
able, requires a more careful analysis, which is performed 
in Sec. [U] For a larger number of dimensions, the one- 
loop bubble exhibits a strong cutoff dependence. 



Appendix D: Calculation of n""''(p,p') 

For our discussion of the role of high momentum (sub- 
horizon) modes to the integral equations (|43| - p5|) in Ap- 
pendix |F] below, we shall need both the mixed (p M 
p') and UV (// < P,p') behaviors of the one-loop bubble 
ti{p, p') for large p'. We analyze the former in detail here, 
specializing to d = 3 spatial dimensions. 

By definition, see Eq. (IB3|) . the mixed contribution 
is such that p < p'. The propagator (jB6|) 

thus reads 



G{p,p') 



f^H,{p)Ht(p'). 



(Dl) 



integral [36l |5 



J^{q,P,p') = J dqG{qp,qp') = 

- qp'H,{qp)Hl_^{qp') - qpHl{qp')H,^,{qp) 



JVPP' 



p2 pr2 



= -Giqp.qp) , 

where we defined the function 

pH^^iip) 



Hu{p) 



(D5) 



(D6) 



The third line of Eq. (jD5|) is particularly useful when it 
comes to approximating the function J-'{q,p,p') for low 
or large momenta. Indeed, the function TZi,(p) has the 
following simple limits: 



7^^(p>l)w^p and 7^^(p<l) 



(D7) 



We are interested in the behavior of the function 
J{p,p') for p <C 1 ^ p'. In the vanishing mass limit, 
V — ?► 3/2, the integral features IR divergences for qp' < 1 
and |1 — g|p' < 1 as can be seen from the low momen- 
tum behavior of the integrand, see e.g., Eqs. (ID23I) and 
(|D24p below. Using a similar technique as in [7, 36], we 
separate the q integral in (ID3p in four contributions^ 



"'0 



1-4 i-oo 



l-JL Jl+ii. 



(D8) 



with /i ~ 1. Reorganizing the various terms and using 
some changes of variables, we arrive at 



By choosing q = |q| and r = |e -I- q| as integration 
variables, the loop integral (|38|) can be written as 



n(p,p') = 



A 



247r2 



J{p,p'), 



(D2) 



with 

J{p,p') 



dqG{qp,qp') / drG{rp,rp') 
"'li-gl 

dqH{q,p,p')[H^ + q,p,p') - ^(|1 - q\,p,p')}. 

(D3) 



In the second line we introduced the (redundant) nota- 
tion 



J{p,p') = J2Mp,p'), 



(D9) 



1=1 



where (omitting the implicit dependence on p and p' for 
brevity) 



Ji = / dqH{q){T{l + q)-Til-q)}, (DIO) 

J2 = ~ rdqF{q){H{l + q)+U{l~q)} (Dll) 
Jo 

grab together all potentially IR dangerous contributions, 
dq{H{q)F{l + q)-n{l + q)T{q)} (D12) 



'H{q,p,p') = G{qp,qp') 



(D4) 



for later use. We exploit the fact that, in d = 3, the r in- 
tegration can be performed explicitly using the indefinite 



^ Note that the integral is U V finite for p 7^ p' , see Appendix [GI 
We can safely send the upper bound to 00. 
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and 



dqH{q)T{l-q). (D13) 



In the following, we adopt a similar notation as the 
one introduced in Appendix iBl for two-point functions to 
distinguish the various regimes in q of the functions H 
and J-. For instance we note ^-[{q^pjp') as 

'HiR{q,p,p') for qp,qp'<fi, (D14) 

^mix(g,p,p') for qp<iJ-<qp\ (D15) 
Hvv{q,p,p') for f^<qp,qp', (D16) 

and similarly for T{q,p^p'). With this notation, we have, 
for the IR contributions, 



Ji« / dqniK{q){:F^W^ + q) - - q)} (D17) 

Jo 

and 



J2« - / rf9^IR((7){H,nix(l+g)+H,„ix(l-'Z)}. (D18) 

Jo 

We separate the integral J3 in two contributions from 
intermediate qp fJ' ^ IP' and high p < qp, qp' momenta: 



T Tinix I tUV 

J3 - -'3 + -'3 ) 



(D19) 



with 

^mix 



3 ~ / d(7{'Hmix('?)^mix(l+'?) -'Hmix(l+q)J^mix(q)} 



and 



(D20) 



J- 



3^^^ / dq{nvv{q)J'vv{l + q)-nvv{l + q)J'vY{q)}- 



Finally 



1-4 



dq'Hmix{q)J'mix{l - q)- 



(D21) 



(D22) 



To evaluate the various integrals above, we use the 
leading IR and mix behaviors 



'Hm{q,p,p') 



Gm.{p,p') 



(D23) 



:Fm{q,P,p') ~ ~qnm{q,p,p') (D24) 



and 



HmUq,P,P) ~ \ 7, ——T , (D25) 



iq 2 



•^mix(g,p,p') ~ Hmix(g,p,p') ( — + ) .(D26) 



The UV contribution (jD21[) has no IR problem and can 
be safely evaluated by setting ly -> 3/2. Using 



/2 e'P /I 

Hi(p)^J , Hi(p) ^ Hi(p) [ - -i 

2 \ TTp i 2 2 ^-f^' 

we get, after some calculations, 



y-vviq) 
J^vviq) - 



2 

g««(p-p') 



i /I 1 
l + _ ^ 

q \p p 



qpp' 



1 1 

p p' 



Hp-p') 

We thus have (recall that qp,qp' ^ 1) 

^uv(g)-7^uv(i + g) - 'Huv(i + q)J^vv{q) 

^g«(2g+l)(p-p') 2gf + 1 + i (^i - 



(D27) 

(D28) 
(D29) 



App'{p' — p) 

jg-i(2g+l)p' 2 



[9(9+1)]^ 



(D30) 



Notice that the high momentum q~^ behavior shows that, 
as expected, the integral J{p,p') is UV finite for p ^ p' . 

Using the various approximations above in the relevant 
integrals and isolating the dominant contributions in the 
limit p <C 1 <C p', we get, after some manipulations. 



J1 + J2 ~ iGmip,p) — 

(p'} JO 

.Gik{p,p) e-'P 



dx 



sm X ~\- X cos X 

^1-26 



{p 



a1 + £ 



(D31) 



and 



J4 « iGm{p,p) 



2p' 

Gmip,p) e'^P' 
e p' 



dq 



[9(1 - 9)1 



1-e 



(D32) 



The contributions J™'^ and J^^ are suppressed by rela- 
tive factors ep'^/p' and ep'^/p''^ respectively. The IR con- 
tribution Ji -|- J2 cancels the l/{p'y contribution from 
J4 and we obtain, finally. 



J{p,p') 



. Gik{p,p) e-'P' 
e p' 



(D33) 



We conclude that the mixed one-loop bubble n™"*^ as- 
sumes the factorized form 

n--(p,p')«^(p)s(p'), (D34) 

where the low momentum piece A{p) is the one character- 
izing the IR contribution (|C8|) with, here, k-I-1/2 ~ 2 — 2e 
and TTp = — A/487r^e, and the high momentum piece is 
given by 



B{p) 



-ip 



(D35) 
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Appendix E: Calculation of n^^(p,p') 

We shall also need the high momentum behavior of the 
one-loop bubble 11^^. We proceed along similar lines as 
in the previous section. As in the case of 11^^ discussed 
in Appendix [C] one needs to regulate the integral in the 
UV since it diverges for -p' — ^ p. We employ the same 
regulator as before, namely a sharp cutoff on physical 
momenta. This essentially'' amounts to cutting the q 
integration with a sharp cutoff A/max(p,p'). 

We shall see that n^^(p,p') has a singular l/|p — p'| 
structure when p' — > p. For our present purposes, this 
function is essentially involved under integrals such as 
(I39p which, we assume, are dominated by this singular 
structure. Therefore, we seek the leading singular contri- 
bution to fiu^(p,p'). 

In the present case, where > 1, potential IR con- 
tributions to the momentum integral in (|D3I) come from 
the regions (taking p < p') q < fi/p' and |1 — < n/p' as 
well as /i/p' < g ^ /x/p and n/p' < |1 — (?| ^ fJ-/p- How- 
ever, these give regular contributions in the limit p' — p. 
The remaining contributions can be separated in three 
pieces: 



UV 



dqHijviq)J^Vvi^ - q), 



(El) 



where x = p — p' and A = A/ max(p,p'). 

As in the case of the field correlator (|B16p . it is in- 
teresting to note that the UV behavior of the function 
\/pp'Il{p,p') assumes a Minkowski-like form with phys- 
ical momenta in place of time. The Minkowski result 
^Minic{t—t' , p) for the one-loop bubble — in the Minkowski 
vacuum — is recovered for p,p' ^ 1 and |lnp/p'| = 
\t-t'\ < 1: 



nMmk(At,p) 



A 



96i7r2 1 \At 



-ip\At\ 



-^^'^l^*! sinpAt 
'jAt\ pAt 



where we used A\p — p'\ « A\t — t'\ 



(E7) 



Appendix F: Influence of subhorizon modes 

In this section, we present a complete solution of Eqs. 
(|43l) ^ (|45p including high momentum modes. We show 
that the result (ISS]) for the p component of the resummed 
function / is not modified, whereas the F component 
(|7ip is only modified by a constant multiplicative factor 
of order unity. 

The key observation for the present analysis is that the 
IR and mix behaviors of the one-loop bubble function 11, 
Eqs. (jCSp and (jD34p . have the factorized expressions 



J, 



T = i dq{UvY{q)TvY{l + q)-HvYil + q)Tuviq)} 



and 



J. 



UV 



dqUvvil + q)J'vY{q), 



(E2) 



(E3) 



U'^{p,p') ^ Aip)A{p') 



(Fl) 



and 



nr(P,p') - Aip)BFip') , Ilf^{p,p') = Aip)Bpip'), 

(F2) 

where we recall the relation of the one-point function 
A{p) to the one used in the main text, see Eq. (p7)) . 



which can be safely evaluated by setting ly — ^ 3/2, i.e., 
using Eqs. (|D28|) - (|D30j) . We obtain, for the leading sin- 



gular behavior when p' ^ p (we also take ^ ^ 1) 



tUV 



tUV 



Ai(p — p') ' 



sin(p — p') 



Ai(p — p') p — p' 



(E4) 
(E5) 



whereas J 2^ ~ j]^^ /p receives a 1/p suppression. 

Considering also the case p > p', our final result for 
the leading singular behavior of the function 11^^ is 



A 



-2iA|x| 



96i7r2 



' , (E6) 



In principle, one should consider the fact that the cutoff also 
applies to r = |e + q|: there should be a term 9{A — qp)9{A — 
qp')6{A — rp)9(A — rp') under the q integral. However, the UV 
divergence arises from very high momenta where r K, q and we 
neglect these subtleties here. 



A[p)^^,l^A(^n^^ 



(F3) 



and where the functions Bp,p{p) are obtained from 
Eq. (imsl) : 



Bf{p) = ReB(p)=-V=^ 



cosp 



Bp{p) = -21mB{p) = ~2y/^ 



smp 



(F4) 
(F5) 



Let us first consider Eq. p4)) for IR modes p,p' <C 1 
and separate the s integral in a low momentum and a high 
momentum piece: ~ Ip ~^ Itj° ^i^h fj, ^ 1. Using the 
notation (|Bl[) - (|B4p . we thus write 



dsif-{p,s)n^^{s,p') 



if{p,p') = tfR(p,p' 



ds/f(p,s)nf(s,p')- 



(F6) 



The first and third terms on the right-hand side corre- 
spond to the IR contribution whereas the second term is 
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the contribution from high momentum modes. The for- 
mer involves the functions H^ , which itself receives con- 
tributions from high momentum modes, as we shall see 
shortly, and the latter involves the mix behaviors IT™" 
and which connect IR and UV modes. Applying a 
similar treatment to Eq. (|45|). we get, for the IR behavior 
iP,p' < 1), 



and, for the "xim" behavior (p' < /i ^ p) , 



(F7) 



nT'{p,p') - fi|!-(p,p') 



dslPrip,s)l';'{s,p') 



(F8) 



Here, the functions ir^^ip,p') = fl'^'''{p' ,p) and fl^^ 
are known whereas the functions and are to 

be determined. Applying again the same technique as 
above to Eq. we get, for the former (which involves 
P,P' ^ A^) 

i';'{p,p')=flf{p,p')+ r dsflf{p,s)i';^{s,p'). (F9) 
Jp 



The solution of this equation is given in Eq. ([65|) . Indeed, 
the integral in Eq. for the p component only involves 
momenta s of the order of the external momenta p and p' 
and thus the IR behavior receives no contribution from 
UV modes. 

For the function 1^™, we use the relation ~ 
—I^^^{p',p) and write (so, below, p ^ p ^ p') 

if'^ip^p') = fi';:^{p,p') + r d.siif{p,s)if''{s,p'), 

Jp 

(FIO) 

where we defined 



^TniP^P) - n--(p,p') + / dsU^^^^p, 

Jfj, 

(Fll) 

These equations have a similar structure as those we 
solved when discussing i]^, i.e., Eqs. (gH) and (gS]), with 
the upper bound replaced by ^. We apply a similar tech- 
nique as that used in Sec. IVl which exploits the fact that 
the function 11™"' factorizes; see Eq. (jF2p . In particu- 
lar, Eq. (jFllI) implies that the function 11™^ can also be 
factorized: 



fl^;^^{p,p') = A{p)B,{p'), 



(F12) 



with 



Bp{p) = Bp{p) - I rfs/UV(p, s)Bp{s). (F13) 



The function is entirely determined by UV physics. 
It solves the following equation: 



i;^''{p,p')^u)!^{p,p')+ / d.n;/v(p,.)/^"V(5,p'), 

(F14) 

which, for |lnp/p'| = \t — t'\ <^ 1 reduces to the equa- 
tion for resumming bubbles in the Minkowski vacuum; 
see Eqs. (jE6P and (jE7|. We note, in particular, that 
this equation does not generate any infrared enhance- 
ment factor. Each new bubble generated by the integral 
equation thus brings a genuine coupling factor A. 

Now, using Eq. (jF9|) . one shows that Eq. (jF10|) can be 
solved as, see Appendix [K[ 



ij''-{p,p') = n;5:^(p,p') + / ds ifip, s)n-^(s,p'), 

Jp 

(F15) 

from which it follows that 



if-{p,p')^A{p)Bp{p'), 



with 



A{p)^A{p)+ I dslf{p,s)A{s) 



(F16) 



(F17) 



We note that Eq. (jF17|) is identical to Eq. ([69|) which we 
solved in Sec. IVl We thus have 



A[p)=^J-^A{\nP^-'^ 



P 



P 



ii+i 



(F18) 



with K — V — e. 

Inserting if'^{p,p') = -if'^ip'.p) = -Bp{p)A{p') in 
(|F7)) . we obtain 

&]}{p,p') = A{p)A{p') \l+ r dsBp{s)BF{s)\ , 



(F19) 

where we used the factorization property ()F1|) as well as 
Eq. (|F17p . Similarly, we get, from (|F8|, 



nir\py) = \ Bf[p) + I dsllf'ip, s)Bp{s) } A{p'). 

(F20) 



We emphasize that the integrals in both (jF19|) and (|F20p 
are finite in d = 3 as can be checked from Eqs. (IE6I) . 
(|F4)) . ((F5|) and (jFlSp . All the ingredients entering the 
right-hand side of Eq. (|F6[) have now been determined. 
Using again Eq. (|F9[) . we finally obtain 



lf-{p,p') = ZAip)Aip'), (F21) 
with the renormalization factor 

Z = 1 + 2 dsBp{s)BF{s) 
Jfj. 

dsds'Bp{s)IlYyis,s')Bp{s'). (F22) 
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It is remarkable that Eq. (IF2ip . which fully takes into 
account UV modes, is essentially identical to the solution 
obtained in Sec. |Vl see Eqs. and ([75]), where UV 

modes were neglected. Superhorizon modes only affect 
the final result through a constant, finite renormalization 
factor Z. This important result shows that there is an 
effective decoupling of UV and IR physics in de Sitter 
space. 

The fact that the statistical two-point function Ip (the 
two-point statistical function of the composite operator 
(fi'^) takes a factorized form, Eq. (jF21[) . in the IR is a re- 
markable result too. As already emphasized in Sec. |Vl 
this is reminiscent of the fact that de Sitter IR fluctua- 
tions can be described by a classical stochastic ensemble 
of growing modes |57H59| . Here, we emphasize another 
interesting result, namely the fact that the spectral com- 
ponent Ip also factorizes in the mix regime, see Eq. (jF16l) . 
This is also rooted in the classical statistical nature of IR 
fluctuations in de Sitter and is to be put in parallel with 
the corresponding factorized form of the spectral function 
of the field ip, Eq. (|BT2)) . 

We end this section with a few remarks concerning 
the renormalization factor (jF22[) . The parametric de- 
pendence of the various factors in (jF22[) on the coupling 
constant A is as follows. From Eqs. (|E6I) and (|F14[) we 
get n'^^ - A and ~ A(l + 0{X)). This implies 



that Bp ^ 
thus have 



Bp{l + 0{\)) 
where 7r„ ^ 



see Eq. (|FT31) . Recalhng that 
A/e, Eqs. (|F4| and we 



Z~l^TTp{\+0{X)). 

As discussed in Sec. IVI[ one has either 



(F23) 



•\/A, in the 

case of a vanishing tree-level square mass, or TTp ^ 1 in 
the case of negative tree-level square mass. In any case, 
for small coupling we can approximate Bp ^ sa Bp p and 
write, using ((F4| and (|F5|) . 



Z 



1 - 2ttp 



l-|-27r„ 



ds 

sin 2/i 
2^ 



sin 2s 



Ci(2/i) 



(F24) 



where we used np 



-7rp/2 in d = 3. 



Appendix G: Renormalization 

The present theory is perturbatively renormalizable in 
D = d + 1 = A. The renormalization of the gap equa- 
tion has been discussed in 33]. The main result is 
that it can be made finite by a standard (Minkowski-like, 
i.e., curvature independent) redefinition of the parame- 
ters TO^, ^ and A. Here, we discuss (for the first time to 
our knowledge in the present context) the renormaliza- 
tion of the four-point function Eq. (|20)) . or (|3T|) . 

As in previous sections, we regularize the theory by 
a simple sharp cutoff on physical momenta, i.e. by us- 
ing the regularized propagator G{p,p') — G{p,p')d(A — 



p)9{A — p'). Such a cutoff on physical momenta actually 
breaks the de Sitter symmetry but is compatible with the 
p- representation, i.e., it does not break the affine sub- 
group of the de Sitter group, which is at the root of the 
p- representation [47]. Moreover, the local Lorentz — and 
thus de Sitter — symmetry is restored after renormaliza- 
tion provided one imposes suitable, i.e., covariant renor- 
malization conditions. An alternative would be to em- 
loy de Sitter invariant regulators, such as dimensional 
S, or Pauli-Villars [gtI] regularizations. 
When discussing renormalization, i.e., the removal of 
the regulator, the decomposition (|40p should be modified 
to take into account a local term 

Il{p,p') = iTTS^ip-p')+np{p,p')~-sign^{p-p')Ilp{p,p'), 

(Gl) 

where tt is a (divergent) constant. This expresses the 
fact that ill{p,p') is a one-loop contribution to the four- 
point function, see (pij) . and thus potentially contains a 
(divergent) contribution with the same structure as the 
tree-level (bare) coupling* ~ i5^(p — p'); see (|37l) . Let us 
now give a general analysis of the divergent contribution. 
We recall the expression with d — 3 spatial dimensions 



n(p,p') 



A r G{qp, qp')G{rp, rp') 
6 Jq qr 



(G2) 



with r = ]e + q]. The divergent contribution 11'^''^ comes 
from very high momentum modes 1 ^ /i^ < qp, qp' < A, 
where is an arbitrary (renormalization) scale. Spe- 
cializing to the case p,p' G C"*" for simplicity,^ we thus 
have^° 



H'^-(p,y) 



A r Glyiqp,qp') 

6 J/i^<|q|<A 9^ 

A '•^ 



487r2 



dq e 



-2iq\p~p'\ 



^ g-2'iA|p-p'| _ ^-2ipr\p~p'\ 



967r2 



i\p-p'\ 



(G3) 
(G4) 
, (G5) 



where A = A/max(p,p') and fir — /i^/ max(p,p'). This 
expression agrees with our previous findings Eqs. (|C18[) 
and (|E6)) . 

To extract the divergent contribution in the limit A — )■ 
oo, we use the well-known result 



sin2A]a;] 



— ?> Trd{x) 



(G6) 



* Equivalently, ili(p,p') is the one-loop contribution to tlie two- 
point correlator of the composite field ip^, whose tree-level con- 
tribution is ~ (5jj(p — p'); see II37I I. 

^ The general result has the same structure with the replacement 
\P-P'\ -(P-P')signc(p-P')- 

Note that it is sufficient to approximate r q in order to get the 
leading (logarithmic) divergent behavior. 
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as well as 



Accordingly, we write 



cos2A|a:| cos2/2^|g;| 

^ 2 In + (G7) 

The latter can be shown as follows. For any test function 
f{x), one has 



cos2A|a;| — cos2//r|a;| 



dx f{x 
^ 2/(0) 



dy 



cosy — cos 

y 



= 2/(0) lim 



Ci(^)~Ci(,) 



2/(0) In ^, 



(G8) 
(G9) 
(GIO) 



where we used /tr/A — fir/-^- We made the change of 
variable y — 2 Ax and took the limit A//ir ^ 1 in the 
second line. The third line is obtained by noticing that, 
since the integrand vanishes at y = 0, one can safely 
replace dy = lim^^o ■ Finally, we used Ci(a;) — 

— dt cost/t = Inx + 7 + 0{x^) at small x. We thus 
get the desired divergent contribution as 



In — 6{p — p') 



487r 



Mr 



(Gil) 



We note that, since the UV divergent term is local, 
one can replace p — p' — > —p{t — t'). One recovers 
the usual Minkowski one-loop coupling divergence for 
y/^tl{p,p') nMi„k(i - t\p): 



Tjdiv 



it-t',p)-^^ln-S{t^t'). (G12) 
487r^ Pr 



This guarantees that the divergence (IGlip can be ab- 
sorbed by the usual Minkowski counterterm. Inciden- 
tally, this ensures that the local Lorentz — hence de 
Sitter — symmetry is properly restored in the usual way 
after the UV cutoff is removed. 

We are now in a position to discuss the renormalization 
of the resummation equation (1391) . For this purpose, it 
is convenient to extract the explicit coupling dependence 
of each term. We define 

U{p,p')^X7T{p,p'), i{p,p') = X-'^{p,p'), (G13) 

such that vr ^ A° and 5 ~ . Equation thus reads 

^{p,p') = \^TT{p,p')-iX [ dsTT{p,s)^{s,p'). (G14) 

Jc 

We have, from the above discussion, 

T^{p,p') = iTTdivS^ip - p') + n^{p,p'), 
where TT-^{p,p') is finite and 



(G15) 



■In- 



A 

Pr 



(G16) 



S(p,p') =z3div%(p-p') + ^^^fep'), (G17) 

where ^f{p,p') is finite, as shown below. Plugging the 
expressions (jGlSp and (jG17l) in Eq. (|G14I) . we get 



1 - A-TTdiv 

and, after some simple manipulations. 



(G18) 



{j),p') = XItt^ {p,p') - iXr / dsTT^ {p,s)'^f {s,p'), 

Jc 

(G19) 

where we introduced the renormahzed couphng 



X 



X 



1 - Att, 



div 



1 



48Tr 



^In- 



(G20) 



This is the standard expression for the renormal- 
ized coupling in the large-iV limit, up to possible 
(renormalization-scheme-dependent) finite parts [131 . It 
shows the trivial character of the present theory, i.e., 
that for any positive bare coupling A, the running cou- 
pling Xr 0+ in the IR, i.e., for /z^/A — 0. Reversely, 
demanding that A^ be finite, Eq. (IG20p expresses how 
the bare coupling A evolves as the cutoff is removed: 
1/A = 1/Ar + TTdiv One get the usual Landau pole in 
the UV, the value A^ of the cutoff where the bare cou- 
pling diverges: Al/pt = exp(487r^/Ar). Finally, we stress 
that the renormalized coupling (|G20I) coincides with that 
needed to absorb coupling subdivergences in the equation 
for self-energy S as expected. 

Using the large momentum behavior of the F and p 
components of the function 7r^(p, s — >■ oo), which can 
be inferred from the previous sections, see Eqs. (|D34|) . 
(jD35p and (jE6p . one can show that the contour integral 
in (jG19p does not lead to UV divergence. This is related 
to the fact that this is originally a time integral, and not 
a momentum loop integral. Thus, we conclude that the 
function 3-^ (p, p') is finite when expressed in terms of the 
renormalized coupling A^ . To make link with the original 
functions 11 and / discussed in the paper, we define the 
finite functions 

n'-(p,p') = A.^/(p,p') , r{p,p') = x-'^f{p,p'), 

(G21) 

in terms of which Eq. (jG19p reads 

r ip, p')^w ip, p')-^J^ ds fr (p, s)r (s, p') • (G22) 



This is to be compared to Eq. (f39| . This shows that the 
discussion presented in the main body of the paper (in- 
cluding appendices) readily applies to the renormalized 
quantities defined here provided one replaces bare quan- 
tities A, n, and I by their renormalized counterparts A^, 
n'', and /'', together with the replacement A — > /i^ in 
momentum loop integrals. 
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To close this section, we note that the bare func- 
tion I{p,p') is not finite, but contains a local divergence 

iXrTTdwS^ip - p'), see Eqs. (|G13|) . (|G17|) . and (|G18|) . 
This is, however, not a problem since / is not directly a 
physical quantity but merely the loop correction to either 
the two-point correlator D, see Eq. ([57]) . or, equivalently, 
the four-point vertex function, see Eq. (|20|) . or (|3T|) . One 
easily checks that these physical quantities arc finite after 
renormalization. Using Eqs. ((Glt| . ([GTtI) . (IGTS)) . ((G20| . 



and (jG22[) . one finds that 



ib(p,p') = 



^5^{p-p')+ir{p.p') (G23) 



is finite and thus so is the four-point vertex function ((3T 
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